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I . IMTRODUCTIOH 


The  problem  of  the  cylindrical  antenna  of  appreciable  diameter  has  received 

considerable  attention  from  a number  of  investigators.  Three  approaches  have 

been  used  in  the  literature.  These  are  (a)  Hallen's  asymptotic  solution  to  the 
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approximate  integral  equation,  and  King's  modification  of  it,  (b)  Storer's 
variational  solution  to  the  approximate  integral  equation,^  and  (c)  Schelkunoff 's 

k 

approximate  modal  solution.  This  report  re-forraulates  the  problem  in  terms  of 
a mathematical  model  for  which  an  exact  solution  is  possible.  The  method  of 
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solution  is  outlined.  A solution  using  Rumsey's  reaction  concept,  which  is  a 
straight  forward  approach  to  the  variational  formulation,  is  also  given. 

Interest  is  restricted  to  steady-state,  harmonic  time-varying  fields.  All 
time-varying  quantities  are  represented  in  terms  of  sinors.^  For  a scalar 
quantity, 

V = -/2  lv|  sin(o3t  + a)  = Im(-/2  V e-^*)  , (l-l) 

where  v is  a time-varying  scalar  and  V = Ivle*^  is  the  sinor  of  v.  For  a vector 
quantity, 

£ = Im(V^  E e-^^)  , (1-2) 

where  6 is  a time -varying  vector  and  E is  the  sinor  of  S.  Eq.  (1-2)  should  be 
interpreted  to  mean  that  the  components  of  £ are  related  to  the  components  of 
E according  to  Eq.  (I"l)- 

The  electromagnetic  field  is  represented  in  terms  of  two  vectors,  the 
electric  intensity  and  the  magnetic  intensity,  the  sinors  of  which  are  E and 
H,  respectively.  The  characteristics  of  a region  are  expressed  in  terms  of 
the  parameters  cr , e,  and  p,  called  the  conductivity,  capacitivity,  and  induc- 
tivity, respectively.  Sources  are  represented  in  terms  of  impressed  electric 


current  density  and  impressed  magnetic  current  density,  the  sinors  of  which 
are  denoted  by  J and  M,  respectively.  The  field  equations  in  terms  of  sinors, 
extended  to  include  impressed  currents,  are 

V X H = (jcae  + o-)E  + J , (1-3) 

- V X E = H + M.  (I"^) 


The  right  hand  side  of  Eq.  (l-3)  is  the  sinor  of  the  total  electric  current 
density,  J^,  and  the  right  hand  side  of  Eq.  (1-4)  is  the  sinor  of  the  total 
magnetic  current  density,  M . In  integral  form,  the  field  equations  are 


(1-5) 

(1-6) 


where  is  the  sinor  of  the  total  electric  current  through  the  surface  S 
encircled  by  the  contour  C,  and  is  the  sinor  of  the  total  magnetic  current. 
The  voltage  is  defined  to  be  the  line  integral  of  the  electric  intensity. 

Thus,  the  relation 


V = 


E • dl  , 

Jc 


(1-7) 


gives  the  sinor  of  a voltage  between  the  "terminals"  of  the  curve  C.  The  value 
of  V is  substantially  Independent  of  all  more  or  less  direct  paths  if  taken  in 
a soTorce-free  region  and  if  the  terminals  are  close  together.  The  conduction 
current  is  defined  to  be  the  surface  integral  of  the  conduction  current  density. 
Thus,  the  relation 


I = 


ds  , 


(1-8) 


gives  the  sinor  of  the  conduction  current.  The  value  of  I in  a "wire"  is 
substantially  independent  of  all  more  or  less  tightly  stretched  surfaces  if 
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the  cross  section  is  small.  The  impedance  Z between  two  wires  close  together 
is  defined  to  be  the  ratio  V/l,  where  V is  the  sinor  of  the  voltage  between 
the  wires  and  I is  the  current  in  the  wires. 

As  a mathematical  aid,  it  is  convenient  to  introduce  auxiliary  functions, 
called  vector  potentials,  in  terms  of  which  the  field  can  be  expressed.  For  a 
linear,  homogeneous  media,  if  M = 0,  one  can  write 

H = V X A , (1-9) 

where  A is  called  a magnetic  vector  potential.  Note  that  A is  not  unique, 
but  is  determined  only  to  within  the  gradient  of  a scalar  since  VxVj^  = 0. 

/V 

One  can  also  use  an  electric  vector  potential  for  regions  in  which  J = 0,  but 
this  is  not  necessary  in  this  report.  From  the  field  equations,  E in  terms  of 
A is  given  by 

(j03e  + o^E  = V x V X A . (I-IO) 

A \mique  solution  for  A in  terms  of  impressed  electric  currents  in  a region 
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homogeneous  everywhere  is  the  potential  Integral  solution.  A solution 
appropriate  to  cylindrical  coordinates  is  given  by  A = it,  where  ijl  is  a solution 

g 

to  the  scalar  wave  equation. 

A powerful  method  of  reducing  electromagnetic  problems  to  simpler  mathe- 
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matical  models  is  the  equivalence  principle.-^  In  the  literature,  the  same 
results  are  usually  obtained  through  the  use  of  Green's  functions,  but  the 
equivalence  principle  is  usually  more  direct.  A concise  statement  of  this 
principle  is  as  follows.  Consider  a siirface  S enclosing  the  sources  of  an 
electromagnetic  field,  ^H.  Now  consider  a mathematical  model  consisting  of 
the  same  surface  S and  with  the  same  matter  (07  e,  n the  same)  external  to  S. 

On  S are  impressed  the  surface  currents 

J = n X Hj 


M = E X n. 


(l-ll) 

(1-12) 


3 ^ 


_4- 


with  the  source  in  the  original  problem  removed.  One  now  has  the  same  field 
external  to  S and  zero  field  internal  to  S.  Since  there  is  no  field  within 
S,  one  can  place  any  material  desired  internal  to  S without  changing  the  field 
external  to  S.  Particularly  useful  choices  of  matter  to  place  within  S are 
(a)  free  space,  (b)  a perfect  electric  conductor,  and  (c)  a perfect  magnetic 
conductor.  For  case  (a),  one  can  use  the  potential  integral  solution.'  For 
case  (b),  it  can  be  shown  that  the  electric  currents  produce  no  field,  giving 
a representation  in  terms  of  only  M.  For  case  (c),  the  magnetic  currents  pro- 
duce no  field,  giving  a representation  in  terms  of  only  J. 

In  this  report,  use  is  made  of  the  theory  of  generalized  Fourier  transforms. 
A transform  pair  is  given  formally  by 
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F(w)  = I F(z)  e 
oo 

F(z)  = ~ J F(w)  e*^^^  dw. 
-00 


(1-13) 

(1-14) 


Eq.  (1-13)  is  called  the  transformation  integral,  and  Bq.  (l-l4)  is  called  the 
inverse  trainsformation  integral.  The  existence  of  F(w)  implies  certain  restric- 
tions on  F(z).  It  can  be  shown  that  the  tremsfonns  of  the  field  vectors  exist 


if  the  sources  can  be  contained  in  a finite  volume 
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A theorem  which  finds 


application  in  this  report  is  Parseval's  theorem  for  transforms.  This  is 


-00 


et 

F(z)  G(z)  dz  = ^ J F(w)  G(-w)  dw. 
-00 


(1-15) 


C /sQ, 

As  defined  by  Rumsey,  the  reaction  between  two  sets  of  soxirces,  J ,M 


Ah  Ah 

and  J ,M  , is 


<a,b>  = 


(E  * i - 1 • M ) dv. 


(l-l6) 


where  E^,h'^  is  the  field  produced  by  the  b sources.  If  the  sources  can  be 


contained  in  a finite  volvune  and  if  the  medium  Is  linear,  the  reciprocity 
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theorem  states  that 

<a,b>  = <b,a>  . (1-17) 

The  self  reaction  of  a source  is  that  given  by  Eq.  (I-I6)  for  a = b.  Note 
that  the  self -reaction  of  a circuit  source  is  VI.  If  a set  of  sources  are 
multiplied  by  some  constemt,  the  reaction  is  increased  proportionately,  that  is 

<Ka,b>  = K<a,b>  . (I-I8) 

Rumsey  presents  a method  of  formulating  problems  to  obtain  an  approximate 
solution  which  is  "best  in  the  physical  sense". ^ The  result  is  essentially 
a stationary  expression,  insensitive  to  small  variations  in  a current 
distribution  about  the  correct  distribution. 

II.  FORMDLATION  OF  THE  PROBLEM 

The  cylindrical  antenna  is  considered  to  consist  of  two  conducting  rods, 
separated  by  a small  gap,  and  excited  by  a generator  connected  across  the  gap. 
This  is  illustrated  in  Fig.  1.  The  generator  is  assumed  to  be  physically 
small  so  that  it  can  be  contained  within  the  cylindrical  gap.  For  this  report, 
the  rods  are  taken  to  be  perfect  conductors. 

The  equivalence  principle  is  now  applied  to  the  cylindrical  surface  just 
enclosing  the  antenna,  including  the  gap  generator.  Choosing  to  back  the 
equivalent  currents  by  a perfect  conductor,  one  has  the  problem  exactly  repre- 
sented by  a single  conducting  rod  with  a loop  of  impressed  magnetic  currents 
encircling  what  was  formerly  the  gap.  This  is  shown  in  Fig.  2.  The  self- 
reaction  seen  by  the  magnetic  currents  will  be  the  same  as  that  seen  by  the 
generator  in  the  original  prbblem. 

Another  representation  of  the  problem  is  in  terms  of  both  the  equivalent 
electric  and  magnetic  currents  over  the  cylinder,  with  both  the  conductor  and 


-7- 


generator  removed.  This  is  called  the  free-space  representation,  and  is  shown 
in  Fig.  5.  The  equivalent  electric  currents  of  the  free-space  representation 
are  equal  to  the  conduction  currents  on  the  rod  of  Fig.  2. 

For  the  mathematical  emalysis,  the  coordinate  system  and  dimensions 
shown  in  Fig.  4 are  used.  The  coordinates  are  the  usual  p,  z circular 
cylindrical  ones.  The  total  length  of  the  antenna  is  2b,  its  diameter  is 
2a,  and  the  length  of  the  gap  is  2c.  Only  antennas  symmetrical  with  respect 
to  the  plane  z = 0 are  considered. 

It  is  noted  that  the  problem  is  symmetrical  about  the  z axis.  One  can 
conclude  therefrom  that  the  impressed  magnetic  currents  of  the  equivalent 
problems  must  flow  only  in  the  J^irection.  The  conduction  currents  along  the 
side  of  the  cylinder  of  Fig.  2,  and  consequently  the  impressed  electric  currents 
of  the  free-space  representation,  must  flow  only  in  the  z-direction.  Over 
the  ends  of  the  cylinder,  the  electric  ctirrents  must  flow  only  in  the  -direc- 
tion. If  the  diameter  of  the  cylinder  is  small,  the  effect  of  the  currents 
on  the  ends  of  the  cylinder  is  known  to  be  small.  There  are  two  types  of 
approximations  that  can  be  made  because  of  the  smallness  of  the  "end  effects". 
One  is  to  formulate  an  approximate  equation  neglecting  the  end  currents,  as 
was  done  by  Hallen.  The  other  is  to  formulate  exact  equations  for  the 

approximate  model  obtained  by  considering  the  conductor  of  Fig.  2 to  be 

15 

hollow.  As  pointed  out  by  Schelkunoff,  only  the  latter  type  of  formulation 
can  be  expected  to  have  an  exact  mathematical  solution. 

Primeury  interest  is  in  antennas  having  a small  gap.  For  this  case,  the 
exact  distribution  of  the  equivalent  magnetic  currents  is  unimportant,  the 
solution  being  essentially  the  same  for  all  well-behaved  distributions.  It 
is  assumed  that  in  the  gap  of  the  original  problem,  and  consequently 
in  the  equivalent  problem,  is  a constant.  Applying  Eq.  (1-6)  to  a path  just 


enclosing  the  magnetic  currents,  and  making  use  of  Eq.  (1-7),  it  can  be  shown 
that 


M^(z)  = ^ , |z|<c.  (Il-l) 

This  is  taken  to  be  the  source  in  the  mathematical  problem.  Of  primary 
interest  is  the  input  impedance  Z,  or  the  input  admittance  Y,  of  the  antenna. 
This  is  the  impedance  or  admittance  seen  by  the  source.  It  can  be  expressed 
in  the  following  ways: 


_ 1 _ V _ <s.s>  _ 


“ <s,s>  " 


(II-2) 

where  <s,s)  represents  the  self -reaction  seen  by  the  source.  Since  the  source 
is  expressed  in  terms  of  magnetic  currents,  the  preferable  form  for  a reaction 
solution  is  the  last  one  on  the  right  side  of  Eq.  (II-2).  If  an  exact 
solution  is  obtained,  all  forms  are  equivalent.  There  is  also  some  interest 
in  obtaining  the  current  distribution  along  the  conductor,  since  it  is 
primarily  this  which  determines  the  radiation  pattern. 


III.  THE  TRANSFOBM  EQUATION 

Considering  the  conductor  of  Fig.  2 to  be  hollow,  an  exact  transform 
equation  can  be  obtained  for  the  problem.  From  symmetry  considerations,  it 
follows  that  the  field  will  be  TM  (transverse  magnetic)  to  the  z direction, 
and  independent  of  General  formulas  for  a symmetric  TM  field  from  a 
cylinder  of  currents  are  derived  in  Appendix  A.  These  involve  the  transforms 
of  the  field  components  in  terms  of  the  transforms  of  the  electric  and  magnetic 
currents.  The  magnetic  current  is  taken  to  be  the  known  source,  given  by 
Eq.  (ll-l).  The  transform  of  M^(z)  is  calculated  to  be 
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Substitutlng  this  Into  the  results  of  Appendix  A,  one  has  for  the  transform 
equation  to  the  problem 


E(w)  = F(w)  + K(w)  J(v)  , 


(III-2) 


where  the  known  f\mctions  are 


^ (aVp^  ) J^(aVp2_w2  ) ^ (lH-3) 

(III-4) 


K(w)  = (p^-w^)  (aVp^-w^  ) (a/p^-w^  ) , 


and  the  unknown  transforms  are 

'b 
-b 


J(w)  3^{z) 


dz. 


E(w)  = ( 


''  “b 

+ 1 ) E^(a,z)  e"'^''^dz. 


(iH'5) 

(ni-6) 
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Eq.  (III-2)  is  an  exact  equation  for  the  mathematical  model,  and  can  be 
expected  to  have  a precise  solution.  It  can  be  shown^^  that  if  p is 
considered  to  be  complex  (corresponding  to  a dissipative  media)  all  terms 
of  Eq.  (III-2)  are  regular  in  the  strip  |lm(w)|  < |lm(p)|.  As  indicated  by 
Eqs.  (III-5)  and  (III-6),  the  inverse  transforms  of  J(w)  and  E(w)  exist 
over  mutually  exclusive  intervals  on  z.  The  equation  is  similar  to  those 
previously  solved  by  a modification  of  the  Wiener-Hopf  method  of  solution, 
except  for  one  important  difference.  The  inverse  transforms  of  the  unknowns 
in  Eq.  (III-2)  do  not  extend  over  semi-infinite  intervals  on  z.  For  the 
benefit  of  anyone  interested  in  attempting  a solution  of  Eq.  (III-2),  the 
following  is  an  outline  of  a procedure  which  shows  some  promise.  Due  to  the 
complexity  of  the  problem,  the  author  has  been  unable  to  obtain  an  exact 
solution. 

Suppose  that  K(w)  can  be  divided  into  the  two  functions 


K(w)  = K3^(w)/iC2(w) 


(III-7) 
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such  that 

00 

p 

E(w)  Kp(w)  dv  = 0 , jzl  ^ b,  (lH-8) 

^<Xi 

r°!- 

J(w)  K,  (w)  dw  = 0 , lz|  > b.  (III-9) 

- OO 

If  and  Kg  can  be  found,  then  an  explicit  expression  for  the  solution  can 
be  written,  obtained  as  follows.  Eq.  (III-2)  can  be  put  into  the  form 

Kg(w)  E(w)  - Kj^(w)  J(w)  = Kg(w)  F(w)  , (III-IO) 


where  the  inverse  transform  of  the  first  term  exists  only  over  the  interval 

|z|>b  and  that  of  the  second  term  exists  only  over  the  interval  |z|<;b. 

Transforming  the  inverse  transform  of  Eq.  (ill -10)  over  the  interval  |z|  < b, 

and  performing  the  z integration  then  gives 

00 

= ,-K^  J '■(=) 

^ - OO 


Since  j(w)  is  a finite  transform,  the  current  distribution  in  Fourier  series 
form  is  given  by  the  inversion  formula^*^ 

cA 

= 2b  y~ , |zl<b  . (III-12) 

n=-oo 


This  can  be  reduced  to  a cosine  series  since  J (z)  is  an  even  function.  A 

z 

solution  such  as  outlined  above  would  be  quite  complicated,  but  might  simplify 

for  small  diameter  antennas  with  narrow  gaps. 

Another  method  of  solution  which  appears  fruitful  at  first  consideration 

is  as  follows.  Suppose  that  one  represents  J (z)  by  a Fourier  series, 

Eq.  (HI-12),  with  the  coefficients  J(njt/b)  to  be  determined.  For  the  true 

solution,  E (a,z)  is  zero  for  z < |b|.  For  any  other  value  of  J (z),  a 
z z 

mean  square  error  can  be  represented  by 

^ E^(a,z,  E*(a,z)  dz  = M (III-15) 


-J(r 
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For  the  true  solution^  M is  zero,  that  is  M is  a minimum.  Thus,  M is  to  be 

reduced  to  a minimum  through  the  choice  of  J(njt/b).  Setting  dM/Sj(nn/b)  = 0 

gives  equations  for  the  determination  of  the  j(njt/b).  For  an  exact  solution, 

there  are  an  infinite  number  of  coefficients  to  determine,  giving  an  infinite 

munber  of  simultaneous  equations.  This  condition  is  vhat  is  referred  to  as 

l4 

"non-final  determination  of  coefficients"  by  Sommerfield,  and  is  of  no 

practical  value.  An  approximate  solution  could  be  obtained  by  taking  a 

finite  number  of  terms  in  the  series,  or  even  taking  some  other  fxinctional 

approximation  for  J (z),  but  the  reaction  approach  presented  in  the  next 
z 

section  appears  to  be  better  suited  for  this  purpose. 

IV.  A REACTION  SOLUTION 

Consider  the  cylindrical  antenna  to  be  represented  in  tenns  of  a 
conductor  encircled  by  a loop  of  magnetic  current,  as  shown  in  Fig.  2. 

The  input  admittance,  from  Eq.  (II-2),  is  given  by 

Y = . (lV-1) 

V 

Since  the  magnetic  currents  are  directly  proportional  to  V,  it  follows  that 
2 

<s,s>~  V , and  Y is  Independent  of  the  input  voltage.  Since  the  form  of 
current  distribution  along  the  conductor  is  known  approximately,  one  can  expect 
a variational  formulation  of  the  problem  to  give  accurate  values  of  Y.  The 
equation  yielding  the  best  answer  can  be  derived  in  a strai^t  forward  manner 
by  following  Rumsey's  reaction  approach.^ 

In  terms  of  the  free-space  representation  of  the  problem.  Fig.  5^  the 
total  reaction  seen  by  the  source  is  the  free-space  self -reaction  of  the 
magnetic  currents,  <m,m>,  plus  the  free-space  reaction  between  the  magnetic 
and  electric  currents,  <m,J)  . That  is, 


<s,s> 


<m,m>  + <m,J>  . 


(IV -2) 
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The  term  <m,m>  involves  only  the  magnetic  currents  and  therefore  depends  only 
on  the  dimensions  of  the  gap.  The  term  <m, j>  involves  both  the  electric  and 
magnetic  currents,  and  thus  depends  on  the  dimensions  of  the  gap  and  on  the 
overall  dimensions  of  the  antenna,  that  is,  on  the  dimensions  of  the  conductor 
of  the  original  problem.  Fig.  1.  Substituting  Eq.  (lV-2)  into  Eq.  (IV-I), 
it  is  seen  that  the  input  admittance  is  the  sum  of  two  terms. 


where 


'if  = 

e c , 


_ .<mj_ml 
~ o > 


Y - 


(Iv-5) 

(lv-4) 

(lV-5) 


The  term  Y will  be  called  the  gap  admittance,  and  the  term  Y will  be  called 
g ” ' c 

the  conductor  admittance.  Since  the  magnetic  current  is  taken  to  be  the  source, 
the  Y can  be  calculated  exactly.  The  problem  is  to  determine  the  best 

s 

approximation  to  Y^.  The  quantity  1/^^,  is  what  most  of  the  previous  investi- 
gators have  called  the  input  impedance,  although  the  necessity  of  including 
the  effect  of  the  gap  haa  been  recognized. 

Along  the  cylinder  coinciding  with  the  original  boundary  of  the  antenna 
(just  inside  of  the  magnetic  current),  the  tangential  component  of  the 
electric  intensity  is  zero.  It  follows  that  the  tangential  component  of  E 
along  the  cylinder  due  to  the  electric  currents  alone  is  just  the  negative 
of  that  due  to  the  magnetic  currents  alone.  Thus,  even  thou^  the  value 
of  J is  not  known,  the  tangential  component  of  E produced  by  J is  known  over 
the  cylinder.  The  reaction  between  the  unknown  currents  J and  any  assumed 
distribution  of  electric  currents  ^ can  therefore  be  calculated.  It  is 
precisely  the  negative  of  the  reaction  between  M and  that  is. 
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<x,d>  » - <x,m>  (IV-6) 

for  an  assumed  electric  current  over  the  antenna  bo\mdary. 

Now  consider  the  question  of  how  to  represent  the  current  along  the 
conductor.  One  STOggestion  would  be  to  express  the  current  in  a Fourier 
series,  and  then  determine  the  coefficients  by  the  reaction  approach.  This 
would  again  involve  an  infinite  set  of  simultaneous  equations  in  an  infinite 
number  of  unknowns, being  of  little  practical  value.  A finite  number  of  terms 
of  the  series  might  be  used  for  £m  approximate  solution,  but  again  a sufficient 
number  of  terms  would  entail  a prohibitive  amount  of  work.  It  is  desired  to 
approximate  the  current  with  as  few  functions  with  unknown  coefficients  as 
possible.  This  can  be  done  by  expressing  the  current  in  terms  of  the  modes 
of  a cylindrical  conductor.  Stratton  shows  that  the  only  symmetric  mode 
external  to  a perfectly  conducting  cylinder  travels  unattenuated  with  the 
speed  of  light. There  will  be  a discontinuity  in  the  current  at  the 
gap,  but  if  the  length  of  the  gap  is  small,  the  discontinuity  can  be  taken 
at  z =0.  From  symmetry  conditions,  the  current  must  be  an  even  function  of 
z.  Thus,  for  a solid  conductor,  neglecting  the  currents  over  the  ends  of 
the  conductor,  the  current  will  be  of  the  form 

J (z)  = . (IV-7) 

z 

This  involves  two  adjustable  constants,  A and  B.  For  a hollow  cylinder, 
in  addition  to  the  external  modes,  there  will  also  be  a set  of  internal 
modes,  corresponding  to  the  symmetric  TM  circular  waveguide  modes.  For  small 
diameter  antennas,  these  modes  sire  all  beyond  cut-off,  and  rapidly  attenuated. 
There  is  also  no  discontinuity  in  the  internal  modes  at  z = 0,  since  the 
boundary  conditions  sure  continuous  in  that  region.  Thus,  for* a hollow  cylinder, 

J (z)  = Ae  + Be ^ C cosh  a z , 

z i n n 

n = l 


(lV-8) 


vhere  the  are  the  propagation  constants  for  the  waveguide  nodes. 

The  internal  modes  are  appreciable  only  in  the  immediate  vicinity  of  the 
ends  of  the  conductor,  but  it  is  necer^'  i*y  to  include  at  least  one  of  them 
in  the  reaction  formulation  of  the  problem. to  satisfy  the  condition  J (+b) 

Z 

= 0.  Then  the  approximation  (lV-8)  again  reduces  to  one  involving  two 
adjustable  constants. 

The  reaction  formulation  for<m,  is  now  derived  as  follows.  The 
ctirrent  distribution  is  represented  by 

J « = AJ^  + BJ^  , (lV-9) 

where  J denotes  an  approximate  current  on  the  cylinder  consisting  of  the 
two  partial  currents  and  J^,  with  adjustable  constants  A and  B.  One 
can  enforce  the  conditions 

<u,a>  = <u,j>  = - <u,m>  , (lV-10) 

<v,a>  = <v,J>  = - <v,m>  , (lV-11) 

that  is,  and  look  the  same  to  ^ as  they  do  to  J.  Bqs.  (lV-10)  and 
(lV-11)  involve  only  free-space  reactions  between  known  sources,  and  can 
thus  be  calculated.  Substituting  for  J from  Eq.  (lV-9)  gives  the  pair  of 
equations 

A<u,u>  + B<u,v>  = - <u,m>  , (lV-12) 

A<v,u>  + B<v,v>  = - <v,m>  . (IV-IJ) 

By  reciprocity,  <u,v>  = <v,u).  Solving  for  A and  B by  Cramer's  rule,  one 
has 


2 ’ 

<u,u><v,v>  - <u,v> 


B = 


<u,v><'u.m>  - <u.u><v,m> 

2 

<u,u><v,v>  - <u,v> 


(lV-14) 


(lV-15) 
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The  reaction  of  interest  is  now 


<m,j>  = <J,m>«=<a,m>  = A<u,m>  + B<v,m> 


= <v.v>Cu.m; 


Cu.Y>Cu.m><y.m>  + <u.u><v.m> 


- <u,u><v,v> 


(lV-16) 


This  is  the  best  approximation  to  <m, using  two  adjustable  constants. 
Substituting  Eq.  (IV-16)  into  Eq.  (IV-5)  gives  the  desired  variational 
expression  for  Y • 


V.  EVALUATION  OF  THE  REACTIONS 

Appendix  A gives  general  expressions  for  the  field  from  a distribution 
of  2-directed  electric  currents  and  ^-directed  magnetic  currents.  Although 
all  currents  are  considered  to  be  on  the  cylinder  p-  the  magnetic  currents 
are  taken  to  be  an  infinitesimal  distance  outside  of  the  magnetic  currents. 
From  Eq.  (I-I6),  the  reaction  between  the  magnetic  currents  and  any  source 


of  interest  is 


<m, x>  = -2rta  J M^(z)  H^(a,z)  dz 


(V-1) 


If  is  the  magnetic  field  from  a distribution  on  the  cylinder,  one 
must  take  the  value  of  just  outside  of  the  cylinder  since  there  is  a 
discontinuity  at  = a.  Making  use  of  Parseval’s  theorem,  Eq.  (I-I5), 
one  can  express  Eq.  (V-1)  as 


<m,x>  = - a 1 M^(-w)  H^(a,w)  dw, 


(V-2) 


where  M and  H are  the  transforms  of  M and  H.  From  Appendix  A,  one  has  for 
the  source 

H^(a,w)  = - ^M^(w)  (aVp^-w^  ) (a*/?^  ) , (V-5) 
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and  for  a source  J^, 


H^(a,w)  ^(w)(p2-w^)  Hp^(aV?^  ) J^(a/  p^-^,2 

First,  consider  the  self -reaction  of  the  magnetic  currents.  The 
general  expression  for  this,  from  Eqs.  (V-2)  and  (V-3),  is 


<m,m)=  J M^(w)  M^(-w) 


) J.  (a-i/^^7  ) dw. 


Substituting  for  M from  Eq.  (III-I),  one  has  for  Y , Eq.  (V-4), 

-A  S 


y - m?-- 

g “ 2 


J"  (aV'^^-w2.)  J^(ai/?^ 


(v-5) 


) dw  (V-6) 


This  is  valid  for  any  size  cylinder.  For  a thin  antenna,  one  can  use  the 
approximate  formulas  for  Bessel  fvmctions  of  small  argument  in  Eq.  (V-6). 


This  gives 


Y »s 
g 


2 


QO  ^ 

Jsin_wcl‘ 

. J 


dw  = JCQ6 


(V-T) 


It  is  noted  that  this  is  just  the  capacitive  susceptance  of  a capacitor 
having  the  dimensions  of  the  gap.  This  result  is  to  be  expected  since 
circuit  concepts  apply  if  the  gap  dimensions  are  small  compared  to  wavelength. 

Now  consider  the  reactions  <m,x>  for  a ciuTrent  along  the  cylinder. 
From  Eqs.  (V-2)  and  (V-4),  one  has  in  general 


<x,m> 


r 

J 


sin  wc 


0^-w^)  (a  ) Jjj(aV?^)  dw. 


(V^) 


Again  this  reduces  to  a simple  expression  if  the  gap  dimensions  are  small. 
If  ?(w)  = 0(w““),  a > 1,  as  w ->  00  , one  can  use  the  approximations  for 


small  argument  Bessel  functions  and  let  sin  wc/wc  ^ 1.  This  gives 
<x,m>  = aV  j*(-w)  dw  = 2naV  J?(0)  . 


(V-9) 


This  is  recognized  as  the  reaction  for  a circuit  source,  VI,  which  is  to 


be  expected  if  the  gap  dimensions  are  small. 
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For  the  reactions  between  electric  currents,  if  the  conductor  is  taken 

to  be  hollow,  one  has  from  Eq.  (I-l6), 

> 

J*(z)  E?^(a,2)  dz. 


<x,y>  = 2ita 


(V-10) 


'Z'  ' z' 

where  J*  and  J?  €ure  any  two  distributions  of  electric  current  along 
z z 

the  cylinder  ^ s a.  Again  using  Parseval's  theorem,  Eq.  (1-15),  one 

M 

<x,y>  = a J E^(a»v)  dw.  (V-ll) 


-Od 

From  Appendix  A,  the  transform  of  the  electric  field  from  a distribu- 
tion of  electric  current  is 


E^(a,v)  = " jj(w)  (p^-w^)  H^^^(aVp^-w^  ) J^(a-/p^-w^  ) . (V-12) 

Thus,  the  general  formula  for  the  reactions  involving  electric  currents 
is 


<x,y> 


- 

2ox 


) dw. 
(V-15) 


As  shown  in  Section  IV,  the  current  along  the  cylinder  should  be  of 

the  form  of  Eq.  (lV-8).  The  functional  form  of  the  current  chosen  for  the 

reaction  approach  must  go  to  zero  at  z s +b,  else  the  self -reaction  as 

given  by  Eq.  (V-IJ)  becomes  Infinite.  This  is  to  be  expected  since  it 

would  require  an  infinite  E to  support  a finite  J at  the  ends  of  a hollow 

z z 

cylinder.  A choice  which  consists  of  the  two  external  modes  plus  the 
dominant  internal  mode  is 


J^(z)  a sin  p(b  - |z|), 

Jj(z)  . cos  s(b  - 1^1)  - 


(V-14) 

(V-15) 


where  a is  the  propagation  constant  for  the 
a is  very  large  and 


mode.  For  small  cylinders, 


1 
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(v-i6) 


(V-17) 


j (O)  SK  COS  pb^ 

z 

for  any  antenna  of  appreciable  length.  Defining  the  integrals 
<x,y>  = - (ata)^  , 

one  has  from  Eq.  (V-13), 

I^u  = 8^  r ) J^(aVp^  ) dw,  (V-18) 

-Otf 

^uv  = ^ ^(w)O^-w^)  H^^^(aVp^-w^  ) J^(al£^  ) dw,  (V-19) 

[j^Cw)]^  0^-w^)  ) dw. 


" 8n3 


(V-20) 


These  simplify  to  some  extent  for  small  diameter  antennas,  but  the  details 
are  complicated  and  outlined  in  Appendix  B.  The  real  parts  of  the  I^y 
Integrals  are  independent  of  a if  the  antenna  is  thin. 

Now  substituting  the  various  results  of  this  section  into  Eqs.  (IV-16) 


and  (IV-5)>  one  has  for  the  conductor  admittance 

I,  ■ v/Ja; 


I cos^pb  + I sinSb  - I sin2pb 
uu  w uv  ^ 


(V-21) 


I - I 

vtu  w uv 


The  complete  expression  for  the  input  admittance  to  the  emtenna  is  given 

by  adding  Y , Eq.  (V-7)»  to  Y . Considering  the  I integrals,  it  is  seen 
6 ^ ^ 

that  for  small  pa,  becomes  very  large  compared  with  and 

Thus,  except  in  the  region  of  small  sinpb  (corresponding  to  the  region  of 

resonance),  the  conductor  admittance  is  given  by 

Y^  « sis_eh  ^ (V-22) 

■‘‘uu 

This  is  precisely  the  result  that  is  obtained  from  the  reaction  approach 
using  one  EMijustable  constant,  and  approximating  the  current  by  Eq.  (V-l4). 
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VI.  DISCUSSION 

The  evaluation  of  the  I in  the  reaction  or  variational  solution 

xy 

to  the  cylindrical  antenna  is  a difficult  task.  As  shown  in  Appendix  B, 
the  real  parts  are  for  all  practical  purposes  independent  of  the  antenna 
diameter,  and  are  fairly  simple  to  evaluate.  However,  the  imaginary 
parts  are  complicated  functions  of  both  Ba  and  pb,  and  no  simple  expressions 
have  been  obtained  for  them.  Evaluation  by  numerical  integration  is  possible, 
but  would  be  an  arduous  job. 

A comparison  of  the  solution  presented  in  this  report  to  previous 

3 

solutions  can  be  made.  Consider  Storer's  variational  solution."^  This 
solution  used  the  same  approximate  integral  equation  as  used  in  King's 
solution.  In  Appendix  C,  it  is  shown  that  Storer's  solution  can  be  put 
into  the  same  form  as  Eq.  (V-21),  except  that  the  I are  slightly  different. 
They  are  given  by  Eqs.  (V-l8)  through  (V-20)  with  J^(a/p^-w^  ) replaced 
by  unity.  This  difference  would  show  up  only  in  the  imaginary  part  of  I^, 
since  the  integrands  of  the  other  two  integrals  become  negligible  before 

differs  appreciably  from  \inity.  It  is  difficult  to  predict  quantitatively 
what  the  difference  in  the  two  solutions  would  be,  but  any  difference 
would  be  evidenced  only  in  the  vicinity  of  resonance.  The  source  of  this 
difference  is  due  to  Storer's  approximation  of  the  Integrand  of  the  starting 
equation,  which  is  inaccurate  in  the  vicinity  of  the  ends  of  the  conductor. 
Thus,  the  "end  effects"  are  obscured,  and  it  is  no  longer  necessary  to 
require  J^(z)  to  be  zero  at  z = +b.  Instead  of  assuming  Eq.  (V-15)  one 

could  now  tEike  jT(z)  = cos  p(b  - |z|),  without  affecting  the  evaluation  of 

z* 

I . Storer  gives  a comparison  of  his  solution  to  that  of  King,  and  shows 
xy 

that  the  two  agree  reasonably  well. 

There  are  some  experimental  values  of  antenna  impedance  in  the  region 
of  resonance  available  which  indicate  that  both  Storer's  and  King's  solu- 
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17 

tion  are  reasonably  accurate.  ' Before  attempting  numerical  calculations 
for  the  solution  of  this  paper,  it  is  hoped  that  a simpler  method  of 
evaluating  the  imaginary  parts  of  can  be  devised. 
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APEEIHDIX  A 


le  Field  from  a Cylinder  of  Currents 


It  was  shown  in  Section  II  that  the  problem  of  the  cylindrical  antezina 
can  be  reduced  to  a cylinder  of  impressed  currents  radiating  into  free- 
space.  It  was  also  concluded  from  symmetry  considerations  that  the  magnetic 
currents  flow  only  in  the  <j>  direction  and  the  electric  currents  only  in  the 
z direction  (if  the  conductor  is  taken  to  be  hollow).  Such  a configuration 
gives  rise  to  a field  that  is  IM  to  the  z direction. 

A circularly  symmetric  TM  field  in  a homogeneous,  source-free  region 
can  be  expressed  in  terms  of  an  A potential  having  only  a z component 


given  by 


»=  ^ J f(v)  dz. 


(A-1) 


The  f\mction  is  a solution  to  Bessel's  equation,  being  if  the  axis 

(2) 

^ = 0 is  included,  and  ' if  the  region  is  included.  The  fianction 

f(w)  is  to  be  determined  from  boundary  conditions. 

Consider  now  a distribution  of  ^-directed  magnetic  currents  and  z- 

directed  electric  currents  on  a given  cylinder  ^ = a.  Let  denote 

the  surface  density  of  the  magnetic  current,  and  J denote  the  surface 

z 

density  of  the  electric  current.  Let  a superscript  e denote  a quantity 
external  to  the  cylinder,  and  i denote  one  internal.  The  boxmdary  con- 


ditions at  the  given  cylinder  are 


- Hj(a,z)  , 
“ E®(a,z)  - E^(a,z)  . 


The  field  components  of  Interest  are 


(A-2) 

(A-3) 


(A-4) 
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2 

= (3  + — ) A^(yo,z)  , (A-5) 

with  A as  given  by  Eq.  (A-1). 
z 

It  can  be  shown  that  the  Fourier  transforms  of  the  field  vectors  exist 
if  the  sources  can  be  contained  in  a finite  volume, or,  more  generally, 
if  the  sources  are  transformable.  Since  the  Fourier  transform  is  unique. 


,w)  = f(v)  (A-6) 

is  the  transform  of  A^(yO,z).  From  Eqs.  (A-4)  and  (A-5)^  it  follows  that 

JcueE^(  ^,w)  = (p^-w^)  A^(^,w)  , (A-8) 


where  E ajid  H / denote  the  transforms  of  E and  H 
z 9 Z I 

(A-3)^  it  follows  that 


From  Eqs.  (A-2)  and 


(A-9) 

(A-10) 


where  J and  M/  are  the  transforms  of  J and  M/. 
z p z p 

Now  the  expressions  for  the  quantities  of  interest  csm  be  written  out 
explicitly  for  the  two  regions.  Denoting  the  unknown  f for  the  internal 


region  by  f^(w),  one  has 


/>,w)  = f^(w)  ) » (A-11) 

2 2 

E^(  /),v)  = ^ f^(w)  . (A-12) 

5j(y0,w)  = f^(w)  J3^(/5i/p^  ) . (A-15) 

Denoting  f for  the  external  region  by  f^(w)  one  has 

^,w)  = f®(w)  p ^ (A-li^) 
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hJ(^,w)  = f®(v)  h[2)(/)V?^  ) • 


(A-15) 

(a-i6) 


One  can  now  determine  the  f fiinctions  from  the  boundary  conditions  at 
jo  B a.  Substituting  into  Eqs.  (A-9)  (A-IO)  gives  a pair  of  equations 

f®(w)  ) - f^(w)  J (al/p^-w^  ) = M,(w)  , (A-l?) 

o S -w  ^ 

f®(w)  ) - f^(w)  J^(aV?^  ) * (A-l8) 

The  determinant  of  the  coefficients  of  the  f functions  is  recognized  as 
j times  the  Wronsklan  of  Bessel’s  equation,  that  is 


) -J  (a-/?^  ) 


(A-19) 


H^2^(aV?^  ) -Jj^Ca  Vp^-w^  ) 

Solving  Eqs.  (A-l6)  and  (A-17)  by  Cramer’s  rule  gives 

f^(v)  = Hi^^\aVp^-w^  ) Mj^(w)  + ) j^(w)j  (A-20) 

'■] 


f®(v)  = ^ 
' ^ 2J 


2 .2  , - , , (a  Vp^-w^  ) j^(w)| 


7^  Jl(ayr-w"  ) M/(v)  + J 


(A-21) 


The  electromagnetic  field  is  nov  determined  at  all  points  in  space. 
Explicitly,  it  is  given  by  substituting  Eqs.  (A-20)  and  (A-21)  into 
Eqs.  (A-ll)  through  (A-16),  and  applying  the  inversion  integral. 


1 


-25- 


APEEMDIX  B 


On  Evaluating  the  I 


Taking  the  transforms  of  the  assumed  currents,  Eqs.  (V-l4)  and 
(V-15),  one  has 


. 2b (cos  wb  - cos 

“-2  2 
P - w 


^(V)« 


2(6  sin  Bb  - w sin  wb)  2(w  sin  wb  + g tanh  gb  cos  -wb] 


2 

p - w 


2^2 
a + w 


2(P  sin  Pb  - w sin  wb)  2(w  sin  vb  + g cos  wbl 


2 2 
P^  - w 


2^2 
a + w 


The  approximation  in  the  second  expression  for  J is  permissible  if  a is 

z 

large,  that  is,  if  the  antenna  is  thin.  The  integrals  of  interest,  from 
Eqs.  (V-18)  throu^  (V-20),  are 

00 

Ixy  = ^ r H^^^(a/p2-w^  ) J^(a  )dw.  (B-3) 

(2) 

All  functions  in  the  integrand  are  real  except  ' . This  can  be 
separated  into  its  real  and  imaginary  paurts  as 

Rejn^^^a/p^-w^)  =J^(a-|/?^)  , lv|  C p , 

ImlH^^^(aV?^)|  = - N^(aVp^  ) , |wl  < p , (B-4) 

= f ^o^*  ^ ^ K1>P  . 

The  Integrand  of  Eq.  (B-3)  is  an  even  function  of  w.  Using  Eqs.  (B-4), 
one  has 

®®^^xy^  = J^(aVp^  ) dw  (B-5) 

" Sp  N^(aV?^  ) J^(a-/p^  ) dw 


2jt^p 


vw 

J (P^-w^)  K^j(aVT^)  Io(a-/^-p^  ) dw 


-26- 


It  is  now  noted  that  for  thin  antennas,  J^(a-/p^-w^  ) ^ 1 for  |w|  < 

Letting  pw  s u,  the  real  parts  of  the  I,  reduce  to  particularly  simple 

xy 

forms.  Noting  that  the  second  term  of  Eq.  (B-2)  is  negligible  compared 


to  the  first  term  for  |w|  < p,  one  has 


Re  (I  ) 
' uv' 


Red  ) 


rl  . 

1 l fcos 
^ 1 J £sin 


cos  tiBb  ~ cos 


1 - u 


- cos 
2 


[sin  6b  - u sin 


1 - u 


- u sin 


1 - u 


(B-T) 


(B-8) 


(B-9) 


These  can  be  readily  evaluated  in  terms  of  sine  suad  cosine  Integrals  or 
by  manerical  integration. 

Unfortunately,  the  imaginary  parts  of  the  have  not  been  reduced 
to  simple  form.  They  can  be  evaluated  using  various  approximations  for 
the  Bessel  functions.  The  resulting  formulas  are  quite  complex,  and  will 


not  be  given  in  this  report. 
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APPENDIX  C 

Alternate  Derlva-^lon  of  the  Reactions 


Instead  of  solving  for  the  field  from  a cylinder  of  currents  as  a 
■boundary  value  problem,  as  vas  done  in  Appendix  A,  one  can  use  the 

7 

potential  integral  formulation.  For  a cylinder  of  z-directed  electric 

currents,  E is  given  by 
z 

“Jpr" 

dj!^'  dz«  J^(z')  , (C-l) 

o -t 

where 


r"  = 1/ a^  + - 2^ a cos(^-^')  + (z-z')^  . (C“2) 

Using  the  approximation  employed  by  Hallen,  King,  and  Storer,^  Bq.  (C-1) 
for  = a.  becomes 


E^(a,z)  =51^0^- 


z') 


j-J3V(z-z')^+a^ 

■/(z-z  ’ )^+a^ 


dz' 


(C-5) 


This  equation  is  accurate  except  in  the  vicinity  of  z *=  -b.  It  thus 

effectively  masks  the  "end  effects"  of  the  antenna.  Making  use  of  the 
11 


known  integral, 

r. 


-JpV(z-z')^+a^  _ 


- Oft 


V(z-2’)^+a^ 


dz  . ) , 

V ® 


(c-4) 


one  can  take  the  transform  of  Eq.  (C-4),  giving 


P(w)  = - 1^  jj(w)  O^-w^)  ) . (C-5) 

It  is  noted  that  this  is  the  same  as  obtained  for  the  exact  formulation, 
Eq.  (7-12),  except  that  the  term  J^(aVp^-w^  ) is  replaced  by  unity. 

This  term  is  approximately  unity  for  small  ^a  and  small  w,  giving  a 
discrepancy  between  Eqs.  (C-5)  anci  (V-12)  only  for  large  values  of  w. 

This  is  to  be  expected  since  the  behavior  of  E^(w)  as  |w|-aroD  is 
determined  by  the  behavior  of  E*(a,z)  as  z-^+b,  and  this  effect  was 
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masked  in  the  derivation  of  Eq.  (C-5)«  Since  ET(w)  goes  to  zero  as 

z 

g-pajwl  w->-+aD,  it  is  no  longer  necessary  to  require  = 0. 

Thus,  instead  of  Eq.  (V-15)>  one  could  take 


J^(z)  = cos  S(h-|z| ) j 


(C-6) 


vith  j'^(z)  as  given  by  Eq.  (V-l4).  The  transform  of  jJ  now  becomes 
z z 

j!(v)  = 8(P  Pb  - . sin  «b)  ^ 

“ 

and  the  transform  of  remains  as  given  by  Eq.  (B-1).  The  integrals 

z 

appearing  in  Eq.  (V-21)  become 


(C-7) 


oo 


(c-8) 


for  the  f emulation  of  this  appendix.  A comparison  of  the  I as  given 

xy 

by  Eq.  (C-8)  with  those  given  by  Eq.  (B-3)  shows  that  only  the  imaginary 
part  of  ml£^t  differ  appreciably  for  the  two  formulations.  A 
quantitative  comparison  cannot  be  made  until  numerical  evaluation  is 
attempted. 


